Fermi's golden rule in a mesoscopic metal ring 



Peter Kopietz and Axel Volker 
Institut fur Theoretische Physik der Universitdt Gottingen, Bunsenstrasse 9, D-37073 Gottingen, 

Germany 
(June 1, 1999) 



Abstract 



We examine the time-dependent non-equilibrium current in a mesoscopic 
metal ring threaded by a static magnetic flux eft that is generated by a 
time-dependent electric field oscillating with frequency to. We show that in 
quadratic order in the field there are three fundamentally different contri- 
butions to the current, (a) A time-independent contribution which can be 
obtained from a thermodynamic derivative, (b) A term increasing linearly in 
time that can be understood in terms of Fermi's golden rule. The derivation 
of this term requires a careful treatment of the infinitesimal imaginary parts 
that are added to the real frequency to when the electric field is adiabatically 
switched on. (c) Finally, there is also a time-dependent current oscillating 
with frequency 2to. We suggest an experiment to test our results. 
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I. INTRODUCTION 



Consider a mesoscopic metal ring threaded by a time-dependent magnetic flux <p(t) that 
has a static component and a part that oscillates with frequency to, 

<f>(t) = + W Bm(wt) . (1) 

By Faraday's law of induction, the oscillating part generates a time-dependent electric field 
directed along of circumference of the ring, E{t) = cos(o;t), with amplitude 

eLE u = 2iru>^ . (2) 

Here L is the circumference of the ring, — e is the charge of the electron, and 0o is the 
flux quantum. We would like to know the induced current around the ring. In the limit 
to — > this is just the usual persistent current PJ0]. But what happens for frequencies in the 
range between 10 8 and 10 13 Hz, which for experimentally relevant rings |3],|4| corresponds to 
Acu< r _1 ? Here A is the average level spacing at the Fermi energy, and r is the elastic 
lifetime. We use units where Ti is set equal to unity. This problem has first been studied by 
Kravtsov and Yudson || (KY), who found that in quadratic order the time-dependent field 
induces (among other terms that oscillate) a time- independent non-equilibrium current Iq 2 \ 
Calculating the disorder average of this current perturbatively, KY found that it has the 
peculiar property that for frequencies exceeding the Thouless energy E c = TiD/L 2 (where 
T> is the diffusion coefficient) the average of Iq 2 ^ does not vanish exponentially, but only 
as uj~ 2 . This is in disagreement with the intuitive expectation that the external frequency 
u leads to a similar exponential suppression of this mesoscopic non-equilibrium current as 
a dephasing rate in the case of the equilibrium persistent current ||[7||. The perturbative 
calculation of KY is based on the assumption of a continuous energy-spectrum, which means 
that the level-broadening due to dephasing, 1/t v , must exceed the average level spacing at 
the Fermi energy, A. If we assume that for low temperature T the dominant dephasing effect 
comes from electron-electron interactions, a simple estimate || shows that 1/t v (uj) < A for 
\u\ < E c in the limit T — > 0. Hence, for frequencies smaller than the Thouless energy the 
spectrum is discrete and the perturbative analysis breaks down. In this work we shall show 
that in this case the term considered by KY is not constant, but grows linearly in time, a 
result which can be understood simply in terms of Fermi's golden rule of time-dependent 
perturbation theory. 

It is important to point out the difference between the current considered here and the 
direct current due to the usual photovoltaic effect. It is well known M that irradiation of a 
medium without an inversion center by an alternating electric field can give rise to a direct 
current (photovoltaic effect). The lack of inversion symmetry can be due to impurities and 
defects in a finite sample. For mesoscopic junctions the photovoltaic direct current has been 
studied in Ref. |10| . In this case the average current vanishes, because disorder averaging 



restores the inversion symmetry. In our case, however, we calculate the direct current 
induced in a mesoscopic ring threaded by a magnetic flux. Because the magnetic flux breaks 
the time-reversal symmetry, the direct current considered here has a finite disorder average. 
Thus, the physical origin of a mesoscopic non-equilibrium current discussed in this work is 
quite different from Ref. 
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II. THE QUADRATIC RESPONSE FUNCTION: WHAT IS WRONG WITH THE 

GREEN'S FUNCTION APPROACH? 



We consider non-interacting disordered electrons of mass m on a mesoscopic metal ring 
threaded by the time-dependent magnetic flux given in Eq. ([!]). Suppose that we have diag- 
onalized the Hamiltonian in the absence of the oscillating flux (i.e. for cf)^ = in Eq. ([]])) for 
the given realization of the disorder. The time-independent part of the Hamiltonian is then 
Ho — J2 a £ acl l c a , where e a are the exact electronic eigen-energies for fixed disorder, which 
are labeled by appropriate quantum numbers a. The operators c£, create electrons in the 
corresponding eigenstates \a). If we now switch on the time-dependent part of the field, the 
Hamiltonian becomes H = H + V(t), with 

fr(t) = ^M*)£<«lA.I0>4<* 

a, (3 

Here 5ip(t) = (<j> u /<j)o) sin(u;i), and P x = —id/dx+ (2tt / L) ((ft / (fi ) is the x-component of 
the one particle momentum operator. As usual, the coordinate along the circumference is 
called the x-direction, and we impose periodic boundary conditions. Using standard non- 
equilibrium Green's function methods, the contribution to the non-equilibrium current that 
is quadratic in the external field is easily obtained [111: 



(mLy 

x K^(oo l} uj 2 )e^ 1+UJ2)t , (4) 

where (p^ is the Fourier transform of the time-dependent part of the flux ([!]) in units of the 
flux quantum (i.e. (f>(t) — <p = 0o / duj'8(p UJ ie~ VjJt ) and the response function K^(ui,U2) is 
given by 



A' :2 V,.-^! X 

a/37 

/(e 7 ) - f{e ) f(e p ) - f(e a ) 



Q(37 e 7 - e a + iQ X + u 2 + iO 



X 



e 1 - Ep + u 2 + iO Ep - e a + ui + iO 



(5) 



with 



P a ^r =<a|^)</3|A > |7)<7|A,|a). (6) 



Here /(e a ) — ( c a c a) is the occupation number, which in a grand-canonical ensemble is the 
Fermi function. Keeping in mind that the time-dependent part of the flux ([]]) corresponds 
to 

<W = TrriW + w) - S(u' - 00)] , (7) 
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it is clear that in this case Eq.(|4]) contains not only oscillating terms, but also a time- 
independent contribution, 



where 



A, 



-e)(27T0 a 



4(Lm) 3 ^ 



3^? ' 



(8) 



(9) 



and 



a/37 £ 1 £a l ® 



/(e 7 ) - f(e p ) f(e ) - f(e a ) 



_e 7 — Sj3 — uj + iO ep — e a + uj + iO 
Defining retarded and advanced Green's functions, 

Gj(e) = - * . n , Gi(E) ' 

£ — E a + ZD 



£ — £ r 



z'O 



(10) 



(11) 



Eq. (|10"D can also be written as 



Pi 



aP"/ 



X 



{poo 
/ def(e + u>)[G2(e + u>)Gf{e)G*(e + u>) 
J oo 



-G A a {E + uo)G a {e)G*{e + co)} 



dE[f{E + U0)-f{E)] 



x [GZ(e + u))G%(e)GA{e + u 
'p( £ )Gj 



-G«(e + u)G*(e)G*(e + u)) 



(12) 



The structure of the Green's functions agrees with the one given by KY in Ref. [|13||. Note 



however, that these authors work in a different gauge: they represent the electric field by a 
scalar potential, so that their expressions contain only a single current vertex. The intro- 
duction of Green's function is useful for calculating disorder averages. It is common wisdom 
that for the calculation of the disorder average of Eq.(|l^) the terms involving products of 
only retarded or only advanced Green's functions can be neglected [|IIJ. In this approxima- 
tion a perturbative calculation of the disorder average of Eq. (|T2|) has been given by KY ||, 



with the result that the associated time-independent part of the non-equilibrium current is 
proportional to u~ 2 for frequencies larger than the Thouless energy. As explained in Sec.|, 
for frequencies u < E c the perturbative expansion is not controlled anymore since the energy 
spectrum becomes discrete. In fact, it will turn out, that the physical behavior is completely 
different in this regime. 
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To demonstrate the break down of the diagrammatic perturbation theory for systems 
with a discrete spectrum, we now show that an exact evaluation of the disorder average 
of Eg. should actually yield an infinite result. Let us therefore go back to the exact 
spectral representation fllCf ) of the response function. Using the formal identity 



p iii5(x) 



x + iO x 

where p denotes the Cauchy principal part, we can rewrite Eq.flTOl) as 

K^(u, -u) = K&(u, -u) + K^icu, -co) , 

with 



a/3 7 £ 1 £c 



e j ~~ e P ~ W 

= -27r 2 ^ReP a/37 [/(e 7 )-/(£ /3 )] 

a/37 

x 5(e.y — e a )8{e 1 — ep — to) . 



The terms with a = 7 in Eqs. (|15|) and (|i~6"D yield the following contributions, 



(13) 



(14) 



(15) 



(16) 



a/3 



d f( £a )-f(s p ) 



f{e a ) - f(ep) 



Ep-uj (e a -ep- coy 



(17) 



K^LJco, -u) = -2n 2 5(0)Y,^P aPa [f(^) - f{ep)\ 



(Xp 



x5(e a -e p -u) . 



(18) 



The right-hand side of Eq.(ITH) is proportional to the infinite factor 5(0). Hence, the term 
Kgg(u, —u) must also be infinite. Because the singular prefactor 5(0) in Eq.(p~8"D does not 
depend on the disorder, this singularity survives disorder averaging ||14j| . Keeping in mind 
that Eq.(|l2D is mathematically equivalent with Eq.(IIO), we conclude that a correct evalu- 
ation of the disorder average (u,—u) must yield an infinite result |pH| . Unfortunately, 
in an approximate evaluation of Eq. QT2D by means of the usual diagrammatic methods this 
5-function singularity is artificially smoothed out, and one obtains a finite result [f|. 



III. ADIABATIC SWITCHING ON 



The infinite term ( Jl8| ) is clearly unphysical. This term is closely related to the in- 
finitesimal imaginary parts iO that have been added to the real frequencies in the spectral 
representation ( |T0| ) for the response function K^ 2 '(u>, —uj). As emphasized by KY |Tj|, the 
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infinitesimal imaginary parts are a consequence of the fact that the response function must 
be causal when the time-dependent part of the Hamiltonian is adiabatically switched on. Let 
us examine the "adiabatic switching on" of the time-dependent perturbation more carefully. 



Following the usual recipe |T5|, we replace the Hamiltonian H + V(t) by H + V n (t), where 
Vri(t) = exp(rjt)V(t). The limit 77 — s- is then taken at the end of the calculation of phys- 
ical quantities. For large enough times t the physical result should be independent of the 
switching on procedure. Indeed, in the appendix we show by explicit calculation that sudden 
switching on produces the same result for the long-time response as adiabatic switching on. 
However, in the latter case one still has to be careful to take the limit 77 — > only after 
the physical quantity of interest has been calculated. We now show that the singularity in 



Eq.(TS) has been artificially generated by taking the limit 77 — > at an intermediate step of 
the calculation. 

By direct expansion of the time evolution operator in the interaction representation 
to second order in the time-dependent perturbation, we obtain the current for adiabatic 
switching on with finite r\ 

(— e)(27r) 2 f°° 
if' it) = ( mL y J dcoid^Scp^Scp^ 

xKS\co 1 ,u 2 )e- i ^ + ^ t , (19) 



with 



K$ (co u co 2 ) = e 2 * £ —. 

^ e 7 - e a + ui + co 2 + 2i7] 



x 



/0 7 ) - /M /M - f( £ a ) 



e 7 - + u 2 + if] Bp - e a + ui + irj 



(20) 



Comparing Eq.(^) with Eq.(|^), we see that the former is multiplied by an extra factor of 
e 2r>t . If we directly take the limit 77 — > 0, this factor is replaced by unity. This is the limiting 
procedure adopted in the usual Green's function approach, where one takes first the limit 
77 — > in Eq.(|20|) and then inserts the resulting expression into Eq. (|19D . In this case we 



recover Eqs.(|8D and (10), which lead to the divergence in Eq.(|18D. We now show that this 
unphysical divergence does not appear if the limit 77 — ^ is taken after the physical current 
has been calculated. Substituting Eq.fl2~0"D into Eq.(|T9"D we obtain 

tf\t) = A„[K${u,, -co) + K$(-u, U ) 

+K%\u,u)e~ 2i » t + K$(-u, -u)e 2 ^] . (21) 

In analogy with Eq. flU]) , we express K^\uj, —00) in terms of products of real and imaginary 
parts 

K$(u>, -co) = K$ p (co, -co) + K$ 5S (co, -co) , (22) 

with 
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K$ p (oj, -u) = 2e 2 * £ ReP a ^[f(e,) - f(ep)\ 

a/37 



X 



Ep-U) 



(e 1 -e a f + {2r ] Y(e 1 -e p -o J y + ^ 



(23) 



2rj 



f{ep)] 



x 



»7 



(£ 7 -e Q ) 2 + (2r / )2(^-^-^) 2 + ^ 2 



(24) 



From Eq.(p3l) it is now obvious that K$ p does not have any contributions from the terms 
a = 7. The finite contribution in Eq.flTTj) is thus an artifact of taking the limit 77 — >- before 



calculating any physical quantities. Let us now focus on the term fl24l) . If we directly take 
the limit 77 — > using 



rj 



hm ,, 



ird(e) 



(25) 



we recover the infinite result (|T8|). However, the structure of the ^-dependent part of Eq.( 
is familiar from the derivation of Fermi's golden rule of elementary quantum mechanics. As 
discussed for example in the classic textbook by Baym [[DJ, terms with this structure should 
be interpreted as a rate, 1.6. ctS Bj contribution to the current that grows linearly in time. It 
is therefore clear that after taking the derivative of Eq.(p4"D with respect to t we obtain a 
finite result if we then let 77 — > 0. A simple calculation yields 



-2 lim y P c 



a p -e p - uy + t 

-2-K^P a p a [f{e ct ) - f{ep)]5(e a -e p 

a/3 



(26) 



Because this expression contains only a single 5-function, after averaging over disorder it 
becomes a smooth function of to. We conclude that to quadratic order in the field the 
non-equilibrium current induced by the time-dependent flux (|I|) has the following three 
contributions, 



I«(^liml«(f) = 4 3) + * 



dL 



(2) 
kin 



dt 



(27) 



where the time-independent part is given by 



4 2) = K lim[^(a;, -u) + K% p (-u, u) 



-(2) 



'(2) 



24- E 



ReP a 



Pj 



X p 



/(g 7 ) - f(£p) 
£~f — Bp — 10 



+ (u) — > —a;) 



(28) 
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The coefficient of the term linear in time is 



= -27rA,£P„„,»[/(e,») -/(£„)] 

x [5(e a -e p ~uj) + {uj ^ -u)] , (29) 

and the oscillating part is 

Jg>(f) = A u lim[^?(o;,a;)e- 2 ^ + K$(-u, -co)e 2 ^] . (30) 

Thus, a time-dependent electric field with frequency uj induces in quadratic order three 
fundamentally different currents, (a) A time- independent contribution 7 t ^ ; as shown in the 
next section, this contribution can be derived from a thermodynamic calculation, (b) A 
contribution tdlj^l/dt which increases linearly in time; this term can be understood in terms 
of the usual golden rule of time-dependent perturbation theory, (c) Finally, there is also a 
time-dependent contribution I@l oscillating with frequency 2uo. When this term is averaged 
over a time- interval larger than cu -1 , its contribution to the current is negligibly small. 

From the above analysis it is clear that the contribution that is proportional to t cannot 
be calculated within the usual Green's function machinery, because in this approach the limit 
?7 — > is taken at an intermediate step of the calculation, causing an unphysical divergence. 
To further support the correctness of the limiting procedure adopted here we show in the 
appendix that Eqs. (|27j-f2l^) can also be re-derived if the perturbation is suddenly (instead of 



adiabatically) switched on. 



IV. THE THERMODYNAMIC ORIGIN OF THE TIME-INDEPENDENT PART 

OF THE CURRENT 

The time- independent part lyf of the non-equilibrium current in Eq.(|2"7|) has been dis- 



cussed by us in Ref. ||16|| . This contribution can be obtained from a thermodynamic calcu- 
lation. In Ref. [16 1 we have assumed (without further justification) the existence of such 



a relation. Let us now put this assumption on a more solid basis. Within the Matsubara 
(imaginary time) formalism one can directly calculate the imaginary frequency version of 
the response function (oji, 0J2) given in Eq.([5]), i.e. 



lUl , IU0 2 J 



X 



e a + iui + 10J2 
- f(e p ) f{e p ) 



f(e a ) 



Bp + iu 2 E(3 — e a + iui 



(31) 



As pointed out by KY [jl3[| , in order to obtain the causal response function, one should 
first continue both frequencies to the real axis with positive imaginary part (iui —>■ u)\ + iO, 
1U2 — > 0J2 ~\~ i0), and then set uj\ = —002 to obtain the constant part of the physical current. 
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On the other hand, if one performs these steps in opposite order (i.e. first sets iuj\ = —iuji 
and then continues itO\ — > uj + iO) one obtains for the current response function 



-UJ) 



Re£ 



a/37 



a(3~r £ ~f £c 



/(e 7 ) - f{e p ) f(e p ) - f(e a ) 



_e 7 — Ep — uj — iO E/3 — e a + u + iO 
Comparing this expression with Eqs.(|l4]-- |T6|) , it is easy to see that 



K^(uj,-uj) = K^(uj,-uj) 



(32) 



(33) 



(2) 

Hence, the time-independent part J th of the current can indeed be obtained from a thermo- 
dynamic calculation JTH]. Note, however, that our analysis of Sec.|TTT| (see also the appendix) 
implies that the terms with a = 7 in Eq.(^) should be omitted from the sum, i.e. the 
physical current is given by 



7(2) 
'th 



A, 



-uj) + kH\-uj,uj) 



(34) 



where 



(2), 
th 1 



UJ, —UJ) 



K£\uj,-uj)-K 



(2) 
diae 



(uj, -uj) 



E 

Ej — Ep — UJ Ep — E a + 



(35) 



see Eq.(fL7|). The direct diagrammatic calculation of the disorder average of ij^ is difficult, 
because the restriction a 7^ 7 in Eq.(^) is not so easy to implement. In Ref. |16[ the following 

limiting procedure was adopted: Instead of directly calculating K^ h (uj, —uj), consider the 



generalization of the imaginary frequency response function (|31f) for electric fields with 
finite wave- vector q, which we denote by K( 2 \iuj, —iuj,q). The limit q — ► is taken after 
the disorder averaged current has been calculated. As shown in Ref. [ 16| , in the diffusive 



regime the function K^(iuj, —iuj,q) is a smooth function of q, so that the limit q — > is 
well defined. The so-defined averaged response function vanishes for frequencies exceeding 

||16|| . On the other hand, perturbative averaging of 
1l7|) shows that this term vanishes as 



the Thouless energy as exp(— y \uj\/2E c 
the contribution from the (unwanted) diagonal term 



for large frequencies. This indicates that the above limiting procedure indeed eliminates 



the contribution of the unphysical diagonal term (|T 
non-equilibrium current. 



to the time-independent part of the 



V. CONCLUSION 

In this work we have shown that a time- dependent flux oscillating with frequency uj 

that pierces the center of a mesoscopic metal ring generates to quadratic order three fun- 

(2) 

damentally different contributions to the current: a constant non-equilibrium current J th , 
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(2) 

a current tdl^l/dt that grows linearly in time, and a current oscillating with frequency 2uo. 
As shown in Ref. [TjJ, the disorder average of the constant term l£ vanishes for frequencies 



exceeding the Thouless energy as exp(— y \u>\/2E c ). The calculation of the disorder average 

of the contribution tdlj^l/dt remains an open problem. A direct perturbative calculation by 
means of the impurity diagram technique is not straightforward, because Eq.(|29|) involves 
three matrix elements but only one energy denominator. Therefore this expression cannot 
be simply written in terms of Green's functions. 

The main result of this work is the prediction of a current tdl^/dt increasing linearly 
with time. From the well-known derivation of Fermi's golden rule |H| it is clear that this 
result is only valid in an intermediate time interval. In particular, the calculation of the 
long-time behavior of the non-equilibrium current requires non-perturbative methods. 

One should keep in mind that our calculation has been performed for non-interacting 
electrons in a random potential, so that our results are valid as long as the spectrum of the 
system is discrete. We have argued in Sec.| that at low enough temperatures this should 
be the case for small external frequencies, |u;| < E c . On the other hand, for frequencies 
exceeding E c the spectrum is effectively continuous. In this regime the conventional Green's 
function methods can be used to calculate the direct current, so that the results of KY || 
should be valid. 

Let us also point out that the linear time-dependence of the current is a consequence of 
the discrete spectrum, and is not related to the adiabatic switching on procedure in Eq. (pD| ) . 
In the appendix we show that sudden switching on yields the same linear time-dependence 
of the current. It seems reasonable to expect that for sufficiently short times the constant 
part of the current is dominant |l(J. We would like to encourage experimentalists to 



measure the non-equilibrium response of mesoscopic metal rings to a time-dependent flux 
in the frequency range 10 8 Hz < uj < 10 13 Hz. 
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APPENDIX: SUDDEN SWITCHING ON 

To confirm that the " switching on procedure" outlined in Sec.|m| yields the correct phys- 
ical results, let us consider a harmonic perturbation that is suddenly turned on at time 
t = 0, 

(pit) = + <f> u Q(t) sin(^) , (Al) 
where 0(t) is the step function. To second order in W the induced current is 

I^\t) = 2A w Re J2 Pafriffo) - /(£«)] 

0/37 
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g2io;t gj(e 7 — e a )t 



(e a - e~f + 2uj)(e a — s@ + u) 

I — e i{e 1 -e a )t 

(e a — e 1 ){e a — ep — u) 



2co 



(e a -ep) 2 -uj 2 



i(Ef3-E a +U))t gi(e 7 — Ea)t 

Bp — 5 7 + UJ 



+ {UJ 



-U 



The diagonal term a = 7 is 



a/3 



sin 2 {uit) 



(e a - ep) 2 - uj' 



sm 



t) + sin 



2(£j3-£a 



'-t) 



(e a - e:/?) 2 - a; 2 



+ 



sm 



t) 



E/3 - E a + UJ 



sm 



't) 



£/3 ~ £a — W 



(A2) 



(A3) 



The terms in the last line can be interpreted in the same way as is done in Fermi's golden 
rule |nj by using the identity 



sm 



m 



Ae 



71 



tS(Ae) for t -> 00 . 



(2) 

It is now easy to see that for large times idiagW yields exactly the same linear in time 
contribution as given in Eq.([29|). The terms with no explicit time dependence in Eq. ( |A2[ ) 
can be identified with 1^ in Eq.(p8[). 



(A4) 
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